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Abstract 

We compute the pressure and entropy of hot QCD in the Umit of large number of fermions, 
A'^f ^ A'^c ~ 1) to next to leading order in A'^f. At this order the calculation can be done exactly, up to 
ambiguities due to the presence of a Landau pole in the theory; the ambiguities are 0{T^ / A.^^^^^^ 
and remain negligible long after the perturbative series (in g^Ni) has broken down. Our results 
can be used to test several proposed resummation schemes for the pressure of full QCD. 
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I. INTRODUCTION 



At temperatures large compared to the intrinsic scale of QCD, T ^ Aqcd, the running 
QCD coupling constant is small, and perturbation theory should be a useful tool for studying 
thermodynamic and dynamic properties of the theory. This is true in principle, at least. 

In practice, there are obstacles to the use of perturbation theory. The first is that hot 
QCD has an intrinsically nonperturbative scale, at wave numbers k ~ g^T, or length scales 
/ ~ Xjg^T. (The running coupling should be evaluated at a scale of order T.) The 
presence of nonperturbative physics at this scale has nothing to do with asymptotic freedom 
or the existence of the scale Aqcd! even at temperatures many orders of magnitude above 
Aqcd; nonperturbative physics sets in at a scale only a little lower than the thermal scale, 
g^T ~ T/ log(T/AQCD)- Some observables are not very sensitive to this infrared scale; but 
any observable will be sensitive to this scale at some order in the coupling expansion. For in- 
stance, the Debye screeninglength carries nonperturbative corrections which are suppressed 
by only one power of g 0,131; photon emissivity is believed to receive nonperturbative cor- 
rections suppressed g^ [3|; and the baryon number violation rate arises entirely from 
nonperturbative physics, and is sensitive at leading order 0]. 

The measurables which show thermal effects with the least sensitivity to the nonpertur- 
bative scale are the standard thermodynamic parameters, such as pressure, entropy density, 
and energy density. For these quantities, perturbation theory can determine all coefficients 
up to 0{g^), and that calculation has been carried out jl,0,0]- (It is also possible to com- 
pute the 0{g^\n{g)) term, and this is in progress [^.) With the coefficients of the pressure 
up to 0{g^) in hand, one can diagnose the quality of the perturbative expansion, by seeing 
how quickly successive partial sums converge. Here we meet the other problem of perturba- 
tion theory at finite temperature; even when the perturbative series exists to some order, it 
shows very poor convergence. This problem is distinct from the presence of nonperturbative 
magnetic physics; for instance it also appears in hot QED, which is free of strongly coupled 
IR behavior 0,3. 

Recently this problem has received renewed attention. Three groups have proposed reor- 
ganizations or partial resummations of the thermal perturbative expansion, which are hoped 
to improve the convergence, and therefore the utility, of perturbation theory. Andersen, 
Braaten, and Strickland, following an earlier proposal by Karsch et al. in scalar field theo- 
ries jl^, have proposed to do so by adding the hard thermal loop (HTL) Lagrangian to the 
action at tree level and subtracting it again with a counterterm treated as formally higher 
order, a technique they call hard thermal loop perturbation theory 0, 0, 0|. Blaizot, 
lancu, and Rebhan have proposed a technique based on 2PI $ derivable methods, together 
with a set of app roximations to render the method tractable, also based on hard thermal 
loops 14 , 3; IB • A similar proposal was recently made by Peshier jl3| • 



All three groups have produced predictions for the QCD pressure, both for pure glue and 
for QCD with fermions. Unfortunately, it is difficult to judge whether the techniques are 
successful, since all we know independently about the QCD pressure are the perturbative 
series, which is only useful at uninterestingly weak coupling, and lattice results from the 
Bielefeld group 3|; which are only at such large coupling that it is not clear whether the 



resummation techniques should still work. The resummation techniques have also been 
applied to scalar field theory, but in that case exact results are only known for O(A^) scalar 
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theory at large A^, which is an almost trivial theory. 

For this reason we think it advantageous to have accurate results for the pressure, in a 
theory as close to QCD as possible. In this paper we present a calculation of the thermal 
pressure in QCD with a large number of fermions, A^f 3> A'^c ~ 1, at next to leading order 
(NLO) in a large A^f expansion and for general g'^N{, which is treated as 0(1). This theory 
contains a great deal of the physics of QCD; for instance, the gauge fields are screened by the 
plasma with the same complicated frequency and momentum dependence as in full QCD. 
The perturbative series for the pressure shows poor behavior, as in full QCD. Of course, 
some of the interesting physics of full QCD is absent at NLO in this expansion. This is 
unfortunate; but on the other hand it is what allows us to solve this theory. For previous 
exact results in 6 dimensional scalar theory, see jl^ (who encountered some similar issues, 
regarding Landau singularities, to those present here). 

We should emphasize now that, since all of the physics of large A^f QCD is also present 
in full QCD, any approximation which shows poor performance in large A^f QCD cannot be 
expected to be valid for full QCD either. The reverse need not be the case; an approxima- 
tion which works at large A^f may stumble over the physics which is left out in the large 
A^f expansion. Therefore, an approximation's performance against our results should be 
considered an optimistic estimate of its performance in full QCD. 

In the next section, we will explain what we mean by large A^f QCD. In Section ITTTl we 
describe how the calculation of the QCD pressure is carried out in this theory. Our results 
and conclusions are presented in Section HVl but we very briefly summarize them here. At 
small effective coupling g'^Ni, the pressure is close to the free theory value, and perturbation 
theory works well. As the effective coupling is increased, P/T^ falls at first, but eventually 
rises; it becomes of order its free theory value about where Landau pole ambiguities become 
uncomfortably large. The entropy shows the same behavior. It is difficult to obtain such 
behavior in an ideal quasiparticle approximation, such as the model of Peshier, Kampfer, 
Pavlenko, and Soff j20|. 

II. LARGE Ni QCD 

By large A^f QCD, we mean QCD with the number of colors A'c taken as fixed and 
0(1), but with the number of fermionic species A^f taken large and the gauge coupling g"^ 
taken small, so that the combination ^^^A^f, which serves as an effective coupling strength, 
remains fixed and 0(1). Therefore, in A^f counting, a factor of g^ is treated as 0(1/A^f). 
This procedure is appropriate because it is the effective coupling which determines the 
convergence of perturbation theory. 

In the following it will be useful to define g'^f^ as 



Here Cp is the quadratic Casimir of the representation containing the fermions and dp 
and c/a are the dimensions of the fermionic and adjoint representations. In SU(A'c) gauge 




(2.1) 
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theory but with adjoint fermions, g^f^ = g'^NiN^. Other than overall factors involving Nc, 
all the nontrivial coupling dependence of the pressure and entropy which we will find will 
depend on g^^. Note that g^^ is renormalization point dependent; at leading order in A^f its 
renormalization point dependence is, 



2 dgls ials) 



Pig ,) - ^ = ^ • (2.2) 

There are no corrections involving higher powers of g'^^; all corrections to this relation are 
suppressed by at least one power of 1/A^f. Therefore, one can solve completely for the scale 
dependence of g^f^, at leading order in 1/A^f; 

1 _ 1 ,H,7,)^ (2.3) 



Clearly the theory has a Landau pole at ALandau ~ fie^'^^^^^s . We will define ALandau so 
that the gauge field propagator diverges at = ALandau! under this definition, ALandau = 
^g5/6g67r /g^g^ There is no hope that some strong coupling dynamics somehow generate 
a strongly coupled UV fixed point, because Eq. (j2.3j) is exact at leading order in 1/A^f. 
Therefore the theory technically does not exist. 

This problem is fatal if ALandau ~ T; but for ALandau ^ T we argue that it is not. The 
Landau pole means that the theory as such becomes ill defined; but it is possible to "rescue" 
the theory by introducing either additional heavy degrees of freedom or high dimension op- 
erators, if they arise at a scale < ALandau- For instance, if the operator F^^D^ D"^ F'^'^ / Af^^^^^^ 
is added to the Lagrangian with coefficient greater than 1, the gauge propagator becomes 
nonsingular at all energies, which is sufficient for our NLO calculation to proceed. This term 
ruins the renormalizability of the theory, but in a way which will not appear in the NLO 
calculation of the pressure. What the Landau pole does mean is, that the exact definition of 
the theory is ambiguous, and quantities like the entropy or the pressure are only completely 
well defined after the physics which resolves the Landau singularity has been specified. The 
size of the ambiguity is, however, suppressed by a positive power of (T/ALandau), which for 
the pressure at NLO turns out to be (T/ALandau)^- So provided that ALandau is kept suitably 
larger than T, the relative size of the ambiguity is tiny and can be ignored.^ As g'^^ is in- 
creased, this new ambiguity remains tiny long after the convergence of perturbation theory 
for the pressure has broken down. Therefore, we will forge ahead and compute the pressure, 
without specifying the UV completion; but we will only work with values of g'^gifJ^) for which 
(ALandau/T") > 35, which provcs abundantly sufficient to keep the relative ambiguity in the 
NLO pressure below 1%. 

The other feature of large A^f gauge theory which we should emphasize is, that the non- 
abelian nature of the theory is not very important. This is because the gauge field self- 
interactions carry powers of g without introducing powers of A'f. In particular the gauge 
self-interactions will have no relevance in the calculation we perform here, and first appear, 
in the pressure and entropy, at NNLO (next to next to leading order) in 1/A^f. Therefore 
our results will apply both to QED and QCD. 



^ Similar considerations, in the context of scalar field theory, appear in Drummond et at, |2l[. 
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FIG. 1: Leading order, 0{Nf) diagram, (A), which is a free fermion bubble, and NLO, 0{N^) 
diagrams, (B), which are a gauge boson loop with any number of fermion bubble self-energy 
insertions. (Coupling counterterm insertions, though not shown, should also be resummed.) They 
can be resummed by the usual Schwinger-Dyson method, to give a loop with a resummed gauge 
field propagator. 

III. COMPUTATION OF THE PRESSURE 

In a relativistic field theory with vanishing chemical potentials, the pressure and free 
energy density are opposites, P = —F; so we compute the pressure by computing the free 
energy, which equals the trace of the sum of all IPI vacuum bubble graphs. This is to 
be done at finite temperature, which means in Euclidean space with time extent (3 = 1/T, 
periodic boundary conditions for bosons, and antiperiodic boundary conditions for fermions. 
To get the thermal pressure, we subtract off the result at zero temperature; the difference 
will be finite, though renormalization is still necessary because the result depends on the 
coupling ^fgjj and this coupling runs. 

At leading order, O(A^f), there is only one diagram, a bare fermionic loop, shown in Fig.^ 
(A). Its contribution to the pressure is^ 

Plo = iVf Tr (^^^ 1^ ^ - ^) (3-1) 
Here the trace is over group and Dirac indices, and the fermionic sum-integral symbol means 




^ We use a [++++] Euclidean and [ - +++] Minkowski metric. Four-vectors are written with upper case 
letters, their space and time components in lower case, so if Q is a 4- vector, then g° is its time component, 
q is the spatial vector and q= \q\ its magnitude. The Minkowski continued frequency is typically denoted 
with a different letter, usually uj. 
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We have written this in 3+1 dimensions, but dimensional regularization is imphed. At 
leading order, since the T = value is subtracted off, this is irrelevant, the result is finite 
and equals 

P^o=4:NfdF . (3.3) 

^ 8 90 ^ ^ 

This is the usual free theory value; ANfdp counts the number of degrees of freedom [the 4 

counts the two spin states each for particle and anti-particle; dp is the number of colors, 

dp = Nc for fundamental representation fermions in SU(A^c)]- 

At NLO, there are an infinite number of diagrams, but the structure is extremely simple; 
all contributions look like a gauge boson propagator with an arbitrary number of fermion 
loops, and counterterm insertions, inserted. No other graphs contribute at 0{N^). This is 
because the only way to add a loop to a graph without changing the power of A^f, is to add 
a fermion self-energy loop (or the required counterterm) to a gauge boson propagator; any 
other way to add a loop introduces a g"^ without another factor of A'f, and ~ 1/A^f. This 
is the reason that large A'f gauge theory is so simple; the large A^f limit is a much stronger 
organizing principle than large A"c. 

This set of diagrams can of course be resummed by the standard Schwinger-Dyson trick. 
Therefore the NLO contribution to the pressure is given by 



NLO 



^Tr (^^In [[G.TiQ) + n,T (Q)] - / In [[G.TiQ) + K:L{Q)] 



(3.4) 



where the trace runs over group and Lorentz indices. The group trace is trivial, it gives dA 
[which is A^^ — 1 for SU(A^c) and 1 for QED]. The bosonic sum-integral symbol means 

t-^f. /I^. (3.5) 

n^jj' is the fermionic plus counterterm contribution to the gauge boson self-energy at finite 
temperature, H^^^^, is the vacuum value, and Gq ^ is the free inverse propagator. Our sign 

convention for 11 is chosen so that our 11 agrees with the usual Minkowski, [H -1 11 . The 

same expression for the pressure can be derived from the Luttinger-Ward relation j22| , with 
a comparable amount of work. 

Both vacuum and thermal self-energies are gauge invariant, since their evaluation doesn't 
involve gauge field propagators. They are also exactly transverse, QpU}^^{ff) = 0, and 
together with rotational invariance this ensures that II^j^ — n(^^^ can be expressed in terms 
of two scalar functions L2c 



nX(Q) - U^:,iQ) = P^'^nT(Q) + Q'^-'U^iQ) . (3.6) 

The exact form of P'^'^ and Q'^'^ are not important for our problem, all that matters to us 
is that they are properly normalized projection operators, = 2, = 1, and P^'^{Q) + 
Q^'^iQ) + Q^Q^ IQ^ = v'^'^ (which is 6'^'^, since right now we are in Euclidean space). For 
the vacuum self-energy, the transverse and longitudinal components are the same, and their 
value is well known; 



w::, = {p^-^ + Q^n n^^c , n.ac = -^^^ I in ^ - ) . (3.7) 



127r2 y""]!' 3 
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Im (Minkowski) 
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- - - - > - - -> - >- - - > - - - - Re (Euclidean) 
<^ — - ^ <^ ^j^j^ 

^^\„^^ y A 

FIG. 2: Points where summation occurs (dots on Euclidean axis) become poles whose residues 
must be summed. The dashed contour sums these poles, and can be continued (arrows) to the 
dotted contour, which evaluates the discontinuity across the imaginary (Minkowski) axis. 

Here /I is the MS renormalization point. 

Some gauge fixing prescription is necessary to evaluate Eq. ()3.4|) . though the resuh is 
gauge invariant since it represents a complete treatment to some order in a systematic 
expansion (this can be explicitly checked). If we choose a gauge where Gg ^ factorizes in the 
same way as 11, the Lorentz trace and logarithm will commute. For instance, in Feynman 
gauge, \G^^Y'^{ff) = rj'^'^Q'^, and (dropping the overall group factor) 

Trln[[Go-Y'(Q) + nfir(Q)] = 2 ln(Q2 + + n.ac) + ln(Q' + + n.ac) + ln(Q') -2 MQ^) . 

(3.8) 

The contribution containing IIt comes from P'^'^; the contribution containing IIl comes from 
Q'^"; the positive ln((5^) contribution comes from Q^Q'^ /Q"^, and the —2\n{Q'^) contribution 
comes from the ghosts, which have no self-energy at this order in A^f. In a general covariant 
gauge the Q^Q'^/Q^ term would be shifted by In(^) with ^ the gauge fixing parameter; but 
this contributes equally to the vacuum and thermal contributions and so drops out.^ The 
next task is to evaluate the thermal self-energy. This has been done by Weldon 23] (in 
Minkowski space, but the continuation is elementary). His result is listed in the Appendix; 
unfortunately it involves an integration which to our knowledge cannot be done in closed 
form. 

Because the self-energies Hl, Ht are only expressed as single integrals and not as closed 
form, analytic expressions, at this point we must either make approximations, or turn to 
semi-numerical methods. We want exact results; otherwise there is no point in our looking 
at this theory. So the calculation will have to become more numerical. The difference 
between a sum-integral and an integral is an inconvenient quantity for numerical analysis. 



^ This is because Tj^qO C = j {dq^ /2tt)C for any C which is independent. The easiest way to see this is 
directly in the time domain, though it also follows easily on rotating to Minkowski frequency. The result 
for Pnlo in Coulomb gauge differs by a factor of h\{q^) from the covariant result, but this also cancels 
between vacuum and thermal for the same reason. 
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especially when issues of UV regulation have not all been addressed. At large values of 
Q^, the argument in the integral in Eq. (|3.4p is smooth, and the difference between a sum- 
integral and an integral is exponentially suppressed, by 0(exp(— Q/T)). So it is only the 
moderate region where we need to handle sum-integration. In this region, this can be done 
by re-expressing the frequency sum as a contour integration and continuing to Minkowski 
frequency. Namely, we note that the function cot(g°/2T) has poles at g° = 2mTT, of residue 
2T. So the Matsubara sum can be replaced by a contour integral, 

J2 HQ' + n(g)) = W J- cot(g72r) hq' + n) , (3.9) 

where the contour runs just above the real axis from — oo to +00, and then just below the 
axis back to —00. This is illustrated in Fig. |21 which also shows that this contour can be 
deformed to pinch the imaginary (Minkowski) axis. The result from upper and lower half 
planes are equal, so the contour integral becomes 

ln(Q2 + n(Q)) = - / — coth(tu/2T) Im In (g^ - + U{q, iu + 0)) . (3.10) 

qO=2nnT ^ 

Note that coth(a;/2T) = 1 + 2/(e^/^ - 1) = 2(nb(u;) + 1/2). The T ^ limit is trivial; 
just leave out the Bose distribution and use the vacuum self-energy. Therefore, Eq. ()3.4|1 
becomes 

y V ^ (^K+i]Im In (g2-cj2+nT+n,ac) - |lm In (g'-cj^+n^ac) J 
+ Pb+^ Im In ilm In 



(3.11) 



J- vac 



Here 11 is always understood to be evaluated at {q.,ioJ + 0). Note that Imfl^ 
(o;^— g^)((yfgg/127r)^(ct;^— g^) is nonzero for timelike momenta; the imaginary parts of the 
logarithms are always at least 0(5'gg) away from vr. There are no infrared problems in 
evaluating the above expression; we discuss ultraviolet issues next. 

The thermal and vacuum contributions fail to cancel for two reasons. First, there is 
the factor of n\,{uj) associated with the thermal contribution, which can be traced to the 
difference between and J dq^ . There are no ultraviolet problems associated with this 
piece, since at large oj it is exponentially suppressed by 0(exp(— cj/T)), while at large g 
but moderate cu, the thermal self-energy has an exponentially small imaginary part, Imll ~ 
exp(— (g— to')/2T). Therefore this piece is exponentially ultraviolet safe. 

Next, there are IIt and IIl. They are only power suppressed at large Q, and can poten- 
tially cause more trouble. To understand this region it is best to go back to the Euclidean 
expression,'^ and expand in small T^/Q^. In this small T^/Q^ limit, the logarithms can be 



We should perform pieces linear in nb in Minkowski space; but terms without nb can be performed either 
in Minkowski space, or in Euclidean space-integrating, not summing, on g^-or using some mixture, as 
convenient. 
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expanded, 



2 In 



Q 



'^+nvac+nT 



+ In 



'^+n.ac+nL _ 2nT+nL 2(nT)^+n^ 
Q2+n.ac QHn^ac 2(g2+n.ac)^ 



(3.12) 



Using Eq. ()A1|) and Eq. ()A7|) from the Appendix, the dominant large G term is 



1 r d^Q 2nT+nL _ in^r^gl 

2 J (27r)4 Q2+n,ac " 45 



/ 



{21XY (g2)3 l-^2^1n[(Q2)/;^2]/i2^2 



d^Q 3(g0)2 - g2 1 



(3.13) 



This term is potentially logarithmically divergent. If one performs an angular average over 
the (Euclidean) direction of Q, at fixed then (3(g°)^ — q^) averages to zero, while 

the rest is a function of only; so any regularization or cutoff procedure which respects 
Euclidean invariance will be UV well behaved. However, if we perform the g° integration 
first, at fixed g, and then perform the d^q integration, then because of the logarithm in 
the denominator, the result does not vanish; instead the g° integral gives a result of order 
((7gg)^T^/g^. Therefore, performing the g° (or Minkowski uj) integration first and the d^q 
integration last holds out the potential for fake logarithmic divergences. 

Any sensible way to UV complete the theory, to deal with the Landau pole, will respect 
Euclidean invariance. Therefore, however we deal with very large values, we should do so 
in a Euclidean invariant way. We choose to apply a large momentum cutoff procedure;^ we 
stop the d^Q integration at = oALandau; for a < 1, varying the value of a between 1/4 and 
1/2 to estimate the irreducible ambiguity. Applying a Euclidean invariant high dimension 
operator cutoff works similarly. 

As we just indicated, it is necessary to cut off the integration at a finite value of Q^, lower 
than A^andau- This also requires some extra care with our Minkowski continuation; we must 
compute the great arcs at complex g° which connect the Minkowski and Euclidean contours. 
Our procedure is this. First, we perform the Minkowski integration up to < ALandauV^ 
and q < gmax ~ 25T (the contribution linear in is carried to larger q). The Minkowski 
integration is done by numerical quadratures, using adaptive mesh refinement. Great care 
must be taken close to the light-cone at large q. To speed things up, the integrals required for 
the self-energies are expressed in terms of a handful of integrals of one variable (either [a;— g] 
or [tu+g]), which are tabulated and spline interpolated so only two integrations nest. Then 
we add the great arc from the Minkowski to the Euclidean frequency axis, and integrate 
over Euclidean g° down to (g°)^ + = g^ax- This accounts for all < g^ax- Then we 
integrate over Euclidean q"^^^ < < aAj^^^^^^, performing first the angular integration 
and then the integration over the magnitude of Q^. Note that (2nT+nL) in Eq. ()3.12p 
vanishes on angular averaging at least to 0{T^^). The lowest order term at large Q^, after 
angular averaging, comes from (H'^+Hl)'^ /(Q'^+H 

vac)^ ~ , SO the irreducible Landau pole 
introduced ambiguity is 0(T*/ALandau)' claimed earlier. 

An alternative procedure is to separate the rib and 1/2 pieces of Eq. ()3.11|) and rotate the 
1/2 piece back to Euclidean space. The n^, piece is then performed in Minkowski space, where 



^ The reader may worry that this wiU introduce gauge fixing dependence. It does not, within the class 
of Lorentz and Coulomb gauges, presumably because the self-energy is gauge invariant; recall that the 
nonabelian nature of the theory is not relevant at the order of interest. 
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FIG. 3: V = PNLo/-PNLO,free and 5'NLo/'S'NLO,free plotted against gig, left, and ALandau/?", right. 
The right plot shows that the features indicating the breakdown of perturbation theory-the pressure 
bottoming out and then rising-occur long before the temperature reaches ALandau- 



it is exponentially UV safe; the 1/2 piece gives Eq. ()H.4j) but with the frequency summation 
replaced by a frequency integration; this can directly be done in Euclidean space, rather 
than by the rather elaborate contour just described.^ The answers agree. 

We compute the pressure at a series of values of g^^^, always at renormalization point 



vre 



-7b 



T ~ 1.76T. 



(3.14) 



with 7^; = 0.577215... Euler's constant. This is the value suggested by dimensional re- 
duction [24]; at this renormalization point, the energy density in an infrared magnetic field 
equals its tree level value, / d^x The result for Pnlo is times a pure function of 

Pnlo = ^dA^T' VigUfl,,^)) . (3.15) 

We have pulled out the free theory value; 2dp^ = 2{N^ — 1) counts the number of degrees of 
freedom in the gauge field, and vr^T^/QO is the usual result for a massless, bosonic field. 

One may also compute the entropy density, which is the temperature derivative of the 
pressure; 



5, 



NLO 



dT 45 V' ' 247r2 dg'^^^ 

This can be done from tabulated results for V by spline interpolation, for instance. 



(3.16) 



This simpler procedure was pointed out to me by Ipp and Rebhan. 
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Perturbation theory at yu.=/u,ppe~',/^pfj,ju,ppe' 



Perturbation theory at fj,= 7tT,27tT,47tT 
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FIG. 4: Exact result for -F^lo, scaled by free value, presented in comparison to the results of 
perturbation theory. Each figure shows three choices for the renormalization point in perturbation 
theory; at left, they are within a power of e of /Udr! right, they are within a factor of 2 of 2itT. 
The left plot also shows the "best" perturbative result, as explained in the text. 

IV. RESULTS AND DISCUSSION 

We present our results for the pressure and entropy, plotted against gl^, in Fig. El The 
figure also shows pressure plotted against ALandau/^- As the coupling is increased, V at 
first falls, as expected (for instance in a quasiparticle picture); then it flattens out and rises. 
The figure shows that the minimum of the pressure occurs when ALandau — 480T. This is a 
scale where the Landau singularity is still far off in the ultraviolet; this minimum is a robust 
prediction of the theory, regardless of the UV completion. 

Note that in an earlier (and published!) version of this paper, the numerical results 
were somewhat different (in particular, higher). This was a result of a coding error in 
evaluating Eq. (j3.11|) : to take the imaginary part of the log, we used the arctangent of the 
imaginary over the real piece. For the longitudinal contribution, we forgot to check whether 
the argument was in the principal range of the arctangent function. Andreas Ipp and Anton 
Rebhan pointed out a discrepancy between the earlier result and their independent numerical 
evaluation of the expressions presented here [i^, and I thank them for correcting me. 

These results are difficult to reconcile with the ideal quasiparticle picture of the plasma 
presented by Peshier, Kampfer, Pavlenko, and Soff [l^l, in which the pressure and entropy 
are described by treating the plasma as a gas of free but massive quasiparticles. In this 
model it is unexpected for to rise as a function of g^^, because this requires that 

the quasiparticle masses get smaller with increasing coupling. Therefore we are doubtful 
that an ideal quasiparticle picture can be a good description of the physics of large iVf QCD. 
Since, as we have argued, large Af is a subset of true QCD, we also don't expect it to be a 
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good description of QCD. 

We also present the predictions of perturbation theory for the pressure, plotted against 
the actual value. The perturbative prediction for Plo agrees with Eq. ()3.3p . The NLO 
prediction, extracted from the results of Zhai and Kastening , is 
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(4.1) 



Note that, unlike in QCD, there is in principle no obstacle to computing the perturbative 
expansion to any order in g'^g. Also there are no appearances of ln((7^g), because the 3-D 
dimensionally reduced effective theory is free at order N^. We have evaluated each partial 
sum from this series, at a number of renormalization points (running the coupling from 
using the exact relation, Eq. ()2.3|) ). in order to compare the perturbative expansion (and its 
renormalization point sensitivity) to the exact result. The comparison is shown in Fig. |3 
The figure also presents a "best" analytic result. This is what one gets by performing the 
dimensional reduction step to the highest order known (two loop, in this case) and then 
solving the dimensionally reduced theory exactly, as advocated recently by Laine [i^. The 
effective 3 dimensional theory can be solved exactly because it is a free theory. The procedure 
turns out to be equivalent to evaluating the perturbative series above at the renormalization 
point where the (s'efr)^^^ term vanishes. As for QCD, the perturbative expansion is valid 
at weak coupling but becomes ill behaved long before the temperature reaches the Landau 
pole (or, in QCD, the phase transition temperature). The "best" analytic result works well 
to surprisingly large coupling but eventually breaks down. 

It is beyond the scope of the current paper to analyze large A^f QCD within each of the 
proposed resummation schemes, to compare to the exact results presented here. We leave 
this task to the authors of those schemes. However, to facilitate the comparison, we present 
our numerical results for the pressure at a number of values of g^f^, in Table UJ 

What are the implications of our results for full QCD? Obviously it does not make 
sense to use large A'f QCD directly to model the behavior of real world QCD; after all, in 
the physically interesting range of temperatures, real world QCD has 3 flavors, which is 
definitely not ^ A^c = 3. This was anyway not the point of studying large A^f QCD. Our 
object was to present a testing ground for other techniques. We expect these techniques to 
be accurate at small g'^^ but, possibly, to break down at some larger value. What does that 
imply for their performance in real QCD? As we have tried to emphasize, success at large 
A'f probably does not ensure success in full QCD; but failure at large A'f almost certainly 
means the technique is not valid. But it would also be nice to have a rough way of equating 
a value of g'^^^ to a comparable value of the QCD coupling as, so if a technique works at 
weak coupling but later breaks down, we can estimate its range of validity in full QCD. 

Probably the most reasonable way to equate g^g with as is to compare values of the Debye 
mass. At least, this makes sense in the picture where screening is the dominant physical 
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TABLE I: NLO pressure, normalized to the free value, as a function of g^^. The numerical error is 
less than 2 in the last place shown, and is dominated by numerical issues; Landau pole ambiguity 
is at most 1 in the last digit in all entries. 

process changing the pressure. The Debye mass for each theory, at leading order, is^ 

„„.(,V,»l).f^^ ™,^(QCD).Hi^,fT^ = ?^P^a.T^ (4.2) 

Equating these, for 3 flavor QCD, gives = ISnas. So ctg = 0.11 (relevant for temperatures 
of order lOOGeV) is comparable to g^g = 6, whereas ctg = -3 (relevant closer to the QCD 
phase transition temperature) is comparable to g^f^ — 17. 

In conclusion, we have presented exact results at next-to-leading order (first nontrivial 



had we chosen /x = MdrC^/^, the large A^f expression would be correct through order (g^g)^. 
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order) for the pressure of large A^f QCD. These results should be useful for testing the 
reliability of resummed perturbative techniques, used for full QCD. 
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APPENDIX A: SELF-ENERGIES 



The self-energies IIt, Hl have been evaluated by Weldon 2^; here we present his results, 
rotated into Euclidean space, and present a few terms of the expansion about small T^/Q^. 

Following Weldon, we define two scalar functions, with normalization chosen to agree 
with his usage. 



2gl^ G{Q) ^ m{Q) - [vac] = 2nT(g) + Ti^Q) 



2gt^ H{Q) = n""(g) 



vac 



1 

Q 



^nL(g) 



so that 



29!, %H{Q) 

Q 



gi^[-:^H{Q) + G{Q) 



(Al) 
(A2) 

(A3) 
(A4) 



Weldon's results for these scalar quantities, rotated into Euclidean space and for general Q, 
are 

?0)2+(2A:-g)2T 



G{Q) 
H{Q) 



27r2 
1 

2^ 



00 



dk 



°° dk 



4fc + — In 



2q (gO)2+(2A;+g)2 
2k I In 



(A5) 



4g 



(gO)2+(2A;+g)2 



q^k , l + 4A;V(gO-ig)2 
1 + 4A;V(gO+^g)2 



(A6) 



Both expressions are pure real, despite the appearance of i in the expression for H{Q). To 
obtain the retarded, Minkowski value, replace g° iuo+Q. 

It is possible to expand each quantity in ^ by series expanding the bracketed 
quantity in small k and then performing the k integration. This is expected to give an 
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asymptotic series, since the k integration goes up to oo but with exponentially decaying 
weight. The next to leading order result, in Euclidean space, is 

^ 45(OT +0(r),(A7) 
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